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A comprehensive first-principles theoretical study of the electronic properties and half-metallic 
nature of finite rectangular graphene nanoribbons is presented. We identify the bisanthrene iso- 
mer of the C2sHi4 molecule to be the smallest polycyclic aromatic hydrocarbon to present a spin 
polarized ground state. Even at this quantum dot level, the spins are predicted to be aligned anti- 
ferromagnetically at the two zigzag edges of the molecule. As a rule of thumb, we find that zigzag 
edges that are at least three consecutive units long, will present spin polarization if the width of the 
system is 1 nm or wider. Room temperature detectability of the magnetic ordering is predicted for 
molecules with zigzag edges 1 nm and longer. For the longer systems studied, spin wave structures 
appear in some high spin multiplicity states. Energy gap oscillations with the length of the zigzag 
edge are observed. The amplitude of these oscillations is found to be smaller than that predicted 
for infinite ribbons. The half-metallic nature of the ribbons under an external in-plane electric field 
is found to be preserved even for finite and extremely short systems. 
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Since their recent successful fabricationji graphene 
nanoribbons (GNRs) have been the focus of exten- 
sive experimental and theoretical efforts. GNRs have 
the same unique hexagonal carbon lattice as carbon 
nanotubes (CNTs), confined to a quasi-one-dimensional 
structure. Hence, they share a variety of interesting 
physical characteristics. Experimental evidence of bal- 
listic electronic transport, large phase coherence lengths, 
and current density sustainability^ accompanied by the- 
oretical predictions and experimental verification of in- 
teresting magnetic properties^'^'^iiiSii^, quasi-relativistic 
behavior fiSiiiiii^iii^ and bandgap engineering capabili- 
tie^i^i^i^ mark GNRs as potential building blocks in 
future nanoelectronic devices. Furthermore, due to their 
planar geometry, standard lithographic techniques may 
be used for the flexible design of a variety of experimen- 
tal devicesi^iiSiii^iii in a controllable and reproducible 
manner. 

Despite the aforementioned similarities, there is a dis- 
tinct difference between CNTs and GNRs. Unlike the 
tubular shaped nanotubes, GNRs, which are long and 
narrow strips cut out of a two dimensional graphene 
sheet, present long and reactive edges prone to localiza- 
tion of electronic states. The importance of these edge 
states was demonstrated for the case of armchair CNTs 
which present a metallic and non-magneti o^^'^^ charac- 
ter in their tubular form, but when unrolled into the 
corresponding zigzag GNRs, they were predicted to be- 
come semiconductingSi2£ with a spin polarized^'^'^i^ii^'^ 
ground state. This ground state is characterized by op- 
posite spin orientations of localized electronic states at 
the two edges of the CNR, which couple through the 
graphene backbone via an antiferromagnetic (AF) ar- 
rangement of spins on adjacent atomic sites. 

In a recent study. Son et al^ have shown that upon the 
application of an electric field, an opposite local gating ef- 
fect of the spin states on the two edges of the ribbon may 
occur. The in-plane field (perpendicular to the periodic 
axis of the ribbon) drives the system into a half-metallic 



state where one spin flavor exhibits metallic behavior, 
while the opposite experiences an increase in the energy 
gap. Chemical doping of the ribbon edges was shown 
to enhance this effect resulting in an efficient and robust 
spin filter devicei^ 

Following these studies on the half-metallic nature 
of one- dimensional zigzag graphene nanoribbons, sev- 
eral papers on the electronic properties)^ magnetic 
properties2^i2ii2^i2& and the half-metallic nature^li^ of 
quasi- zero- dimensional graphene-based structures have 
appeared. These concluded that the spin polarized 
character of the zigzag graphene edges persists also in 
nanometer scale islands of graphene, and that the half- 
metallic nature of the structures studied by Son et al^ 
may be an artifact of the level of theory they appliedi^ 
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FIG. 1: The three sets of finite GNRs studied. The notation 
corresponds to the number of hydrogen atom passivating the 
edges. A N X M ribbon has A'' hydrogens on the armchair 
edge and M hydrogens on the zigzag edge. 
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It is the purpose of our paper to present a comprehen- 
sive and systematic analysis of the electronic properties of 
finite-sized graphene nanoribbons. We present a general 
rule of thumb for the existence of a spin-polarized ground 
state in aromatic carbon based materials and identify the 
molecules C28H14 (Phenanthro[l,10,9,8-opqra]perylene) 
and CsgHie (Tetrabenzo[bc,ef,kl,no]coronene) as the 
smallest hydrocarbon structures to possess a magnetic 
ground state. The stability of the magnetic ordering, the 
energy gap dependence on the dimensions of the system, 
and the effect of an externally applied electric field are 
also studied. Unlike previous studies)^ we show that 
the half-metallic nature of the ribbons (whether periodic 
or finite) is robust and insensitive to the level of theory 
used. 

To this end we study a large set of finite rectangu- 
lar nanoribbons of different widths and lengths using 
three levels of density functional approximations. Wc 
notate the ribbons according to the number of hydro- 
gen atoms passivating the edges, such that a, N x M 
finite ribbon has N hydrogen atoms on its armchair edge 
and M atoms on its zigzag edge (see Fig. [1]). The rib- 
bons we consider constitute three sub-sets correspond- 
ing to three ribbon widths: 4 x M, 6 x M, and 8 x M, 
where M = 2,. ..,16. All the calculations presented in 
this work were carried out using the development ver- 
sion of the Gaussian suite of programsi^ Spin polarized 
ground state calculations were performed using the lo- 
cal density approximation^^ (LDA), the semi- local gra- 
dient corrected functional of Perdew, Burke and Ernz- 
erho^i^i^ (PBE), and the screened exchange hybrid 
density functional due to Heyd, Scuseria and Ernzerhof 
(HSE06)j^i^i^i^ which has been tested in a wide va- 
riety of materials and has been shown to accurately re- 
produce experimental band gaps^i^i and first and sec- 
ond optical excitation energies in metallic and semicon- 
ducting single walled CNTsi^i^S The inclusion of short- 
range exact-exchange in the HSE06 functional makes it 
suitable to treat electronic localization effects^i^i^i^i^ 
which are known to be important in this type of materi- 
als 4>^>^i^i^i^i^i^'^>^>^i^>^i^>^i^i^<^>^ This is further 
supported by the good agreement which was recently ob- 
tained between predicted bandgap&i^ of narrow nanorib- 
bons and measured values. To obtain a reliable ordering 
of the different magnetization states we find it important 
to relax the geometry of the finite GNRs for each spin 
polarization. Therefore, unless otherwise stated, all re- 
ported electronic properties are given for fully optimized 
structures for each approximate functional using the po- 
larized 6-31G** Gaussian basis setiH It should be noted 
that since our calculations are performed within a single 
determinantal framework, we can determine only the to- 
tal spin vector projection along a given axis rris and not 
the total spin. 

Given that most organic molecules are known to 
be diamagnetic, while periodio^i^i^i^iii^i^iSi and finite 
graphene clustera2^i2ii^i21 are predicted to have a spin 
polarized ground state, we begin by studying the emer- 



gence of magnetic ordering in ultra-short and ultra- 
narrow rectangular graphene nanoribbons. In panels 
(a)-(e) of Fig. [2] we present the ground state spin den- 
sity of representative structures of each of the subsets 
studied as obtained using the IISE06 density functional. 
All the structures studied have compensated lattices in 
the sense that the number of carbon atoms belonging to 
each graphene sub-lattice are balanced. Therefore, ac- 
cording to Lieb's theorem for bipartite lattices^i^ they 
have no total spin moment. The two smallest struc- 
tures that present magnetic ordering in their ground state 
are found to be the Phenanthro[l,10,9,8-opqra]perylene 
(bisanthrene) isomer of the C28II14 molecule and the 
Tetrabenzo[bc,ef,kl,no]coronene isomer of the CaeHig 
molecule. The later was recently identified as spin po- 
larized using gradient corrected and hybrid density func- 
tionals.— This is a remarkable finding regarding the ex- 
pected aromatic hydrocarbon nature of these two iso- 
mers especially in light of a previous study that predicted 
the ground state to be diamagnetic and the first excited 
state to be of diradical nature (see panel (f) of Fig. (2]).^^ 
Hence, our HSE06 results suggest that the energy gain 
from the antiferromagnetic ordering of spins obtained in 
the hexaradical Clar structure (see panel (g) of Fig. [2]) is 
higher than the aromatic stabilization even for this small 
organic molecule. 




FIG. 2: Isosurface spin densities of the antiferromagnetic 
ground state of the 04x03 (panel (a)), 06x03 (panel (b)), 
08x03 (panel (c)), 06x08 (panel (d)), and 08x07 (panel (e)) 
GNRs as obtained using the HSE06 functional and the 6- 
31G** basis set. The sorting of the atoms into zigzag, arm- 
chair and seam regions is indicated in panel (d). Panels (f) 
and (g) represent the diradical and hexaradical Clar struc- 
tures of bisanthrene 

From a qualitative analysis of the spin density maps we 
conclude that the edge atoms can be divided into three 
subgroups: atoms which distinctly belong to a zigzag 
edge, atoms that distinctly belong to an armchair edge, 
and atoms belonging to the seam region between the 
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two types of edges (see panel (d) of Figl2]). As a rule 
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FIG. 3: Spin density isosurfaces of the 04 x 16 GNR high spin 
multiplicity states as obtained at the PBE/6-31G** level of 
theory. Panel (a); ferromagnetic arrangement of the nis — 3 
state. Panel (b): mixed configuration of the rUs — 1 state 
presenting a short spin wave on one of the zigzag edges. 

of thumb we find that for rectangular ribbons of width 
N = 4 and above, every carbon atom that distinctly be- 
longs to a zigzag edge will be spin polarized. Seam region 
edge atoms will have a lower spin polarization, while edge 
atoms on a pure armchair edge are considerably less po- 
larized. Therefore, the N x 2 structures, which have no 
distinct zigzag edge atoms show a diamagnetic ground 
state, while the iV x 3 (and longer) structures have mag- 
netic ordering in their ground state. 

An interesting finding is the fact that the maximum 
MuUiken spin polarization on the zigzag edge carbon 
atoms seems to be approximately constant for all the 
systems studied and depends mostly on the density func- 
tional chosen. For the local density approximation this 
value is usually between 0.26 and 0.30, whereas the max- 
imum MuUiken spin density is within 0.30 — 0.37 at the 
PBE level of theory. For the HSE functional, the range is 
0.40 — 0.47. A second important conclusion is that for the 
short systems (3 < M < 8) the energetic ordering of the 
spin states is E^i < E-^i < Eq. Here E-^ is the total en- 
ergy of the antiferromagnetic ground state (Fig. [21), E^^ 
is the energy of the ferromagnetic arrangement with a to- 
tal spin moment projection of rris = 1 and parallel spins 
on both zigzag edges (see panel (a) of Fig. [3]), and Eq is 
the total energy of the diamagnetic state. For longer rib- 
bons, while the ground state remains antiferromagnetic, 
the above lying state is not necessarily the rris = 1 state 
anymore. The number unpaired electrons is not sufficient 
to sustain the desired maximum spin polarization on the 
zigzag edges while maintaining the ferromagnetic order- 
ing. Instead, a mixed configuration is obtained where 
both spin polarizations appear on the same ribbon edge 
forming a type of a short spin wave (see panel (b) of 
Fig. [3]). As a result the energy of the — 1 state is 
raised, usually above the nis = 2 state, which in turn 
becomes the lowest energy state with finite total spin. 

To quantify these findings we study the stability of 
the antiferromagnetically ordered ground state with re- 
spect to the above lying higher spin multiplicity state. 
In Fig. 3] the energy differences between the antiferro- 
magnetic ground state and the first higher spin multi- 
plicity state are presented for the three sets of ribbons 
studied. All diagrams are characterized by sharp max- 
ima structures that correspond to the ribbon length at 
which a transition in the magnetization nature of the first 
higher spin state occurs, as discussed above. We note 
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FIG. 4: Energy differences between the antiferromagnetic 
(AF) ground state and the above lying higher spin multi- 
phcity ferromagnetic (FM) or mixed state for the three sets 
of ribbons studied as calculated by the local spin density ap- 
proximation (upper left panel), PBE functional (upper right 
panel), and the HSE06 functional (lower left panel). The 
green line represents fesT at room temperature. Notice the 
different energy scales obtained for the different functionals 
calculations. 



that even though the antiferromagnetic state of 04x03 
(bisanthrene), as calculated at the HSE06 level of the- 
ory, is only 0.003eV (30K) below the above lying ferro- 
magnetic state, the 04x04 structure has an energy differ- 
ence of ~ O.lSeV (1700K), suggesting the detectability 
of its magnetic ordering even at room temperature. Fur- 
thermore, the HSE06 results indicate considerably large 
energy differences between the antiferromagnetic ground 
state and higher spin states for all ribbons with > 4 
and length exceeding Inm. 

Before discussing the effect of an electric field on the 
electronic properties of finite nanoribbons, it is essential 
to study their ground state characteristics in the absence 
of external perturbations. The geometry dependence of 
the HOMO-LUMO (highest occupied molecular orbital 
and lowest unoccupied molecular orbital, respectively) 
gap would be the most important parameter to address. 
In Fig.[5]the energy gap as a function of the length of the 
ribbon are presented, for the three subsets of ribbons con- 
sidered and the three density functional approximations 
used. Even though up to now we have regarded the stud- 
ied structures as zigzag nanoribbons, one can also think 
of them as wide and short armchair ribbons due to their 
finite length. These type of ribbons arc known to present 
remarkable bandgap oscillations as a function of the rib- 
bon widtbi^ii^ for infinitely long armchair GNRs. Such 
oscillations can be clearly seen in Fig. [5] especially for 
the PBE results (upper right panel). The periodicity of 
the oscillation appears to be somewhat different than the 
threefold period obtained for the infinitely long counter- 
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FIG. 5: HOMO-LUMO gap values for the three sets of ribbons 
studied as calculated by the local spin density approximation 
(upper left panel), PBE functional (upper right panel), and 
the HSE06 functional (lower left panel). Notice the different 
energy scales obtained for the different functionals calcula- 
tions. 



parts. This, however, is a result of the fact that in order 
to prevent dangling carbon bonds, the width step taken 
for the finite systems is twice that taken in the infinitely 
long armchair ribbon calculations. The amplitude of the 
oscillations is found to be considerably damped due to 
the finite size of the ribbons in agreement with a previ- 
ous study.^^ As expected, when the length of the arm- 
chair edge {N) is increased, the oscillations amplitude 
increases as well. It is interesting to note that, in gen- 
eral, the HOMO-LUMO gap is inversely proportional to 
the width (N) and the length (M) of the finite GNR in 
accordance with the semi-metallic graphene sheet limit. 
Therefore, in order to obtain energy gap tailoring capa- 
bility one will have to consider GNRs with long armchair 
edges (large N values) and short zigzag edges (small M 
values). This will increase the amplitude of the energy 
gap oscillations while maintaining overall higher gap val- 
ues. 

Having explored the electronic properties of unper- 
turbed finite GNRs, we now turn to discuss the answer 
to the main question addressed in this study, namely, 
will finite nanoribbons turn half-metallic under the influ- 
ence of an external electric field? In a previous study, it 
was argued that the inclusion of Hartree-Fock exchange 
in the approximate density functional results in the dis- 
appearance of the half metallic statci^ This was shown 
to be not true in periodic systems where screened ex- 
change (and unpublished full exact exchange) calcula- 
tions have shown half-metallic behavior similar to that 
obtained with semi-local and gradient corrected function- 
als.-^ In Fig. [S] we present the spin-polarized HOMO- 
LUMO gap dependence on an external in-plane electric 



FIG. 6: Spin-polarized HOMO-LUMO gap dependence on 
the strength of an external in-plane electric field for three 
representative finite nanoribbons as calculated by the local 
spin density approximation (upper left panel), PBE func- 
tional (upper right panel), and the HSE06 functional (lower 
left panel). Fixed geometries of the relaxed structures in the 
absence of the external field at each level of theory were used. 
Notice the different energy scales obtained for the different 
functionals calculations. 



field applied perpendicular to the zigzag edge, for three 
representative finite GNRs. Similar to previous calcula- 
tions for periodic systemSf^i^ in the absence of an exter- 
nal field the a and /3 gaps are degenerate. Upon the 
application of the field, electrons having one spin fia- 
vor experience a smooth increase in the HOMO-LUMO 
gap while the opposite spin flavor experiences a decrease 
in the gap. This gap splitting continues up to a point 
where the decreased gap vanishes creating a degenerate 
zero energy state. As expected, wider systems require a 
lower onset field to obtain this half-metallic statej^ At 
this point, due to the increased mobility of the metallic 
electrons, further increase in the external field results in 
spin transfer between both edges thus reducing the total 
spin polarization and the energy gap splitting. At higher 
electric fields the systems become diamagnetic— All 
three representative structures present the same features 
described above. The main difference between the LSDA, 
PBE and HSE06 results is the zero-field HOMO-LUMO 
gap, and thus the onset electric field required to induce 
half-metallic behavior. Nevertheless, they all predict the 
appearance of the half-metallic state at an appropriate 
electric field strength. 

To summarize, we have presented a detailed study of 
the electronic properties of finite rectangular graphene 
nanoribbons. Bisanthrene is predicted to be the smallest 
organic hydrocarbon structure to present a spin polar- 
ized ground state. This state is characterized with an- 
tiferromagnetic ordering of spins at the to zigzag edges 
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of the molecule. As a rule of thumb we find that zigzag 
edges, that are at least three consecutive imits long, will 
present spin polarization if the width of the system is 
Inm or wider. For systems with zigzag edges 1 nm and 
longer, our HSE06 results predict the antiferromagnetic 
ordering to be considerably stable with respect to higher 
spin multiplicity states, suggesting the room tempera- 
ture detectability of their spin polarization. Depending 
on the dimensions of the system higher spin states can 
have antiferromagnetic ordering or be of mixed nature 
with spin wave characteristics. Similar to periodic rib- 
bons, HOMO-LUMO gap oscillations with the length of 
the zigzag edge are observed. The amplitude of these 



oscillations however, is lower than that predicted for in- 
finite ribbons. The half-metallic nature of the ribbons 
under an external in-plane electric field is preserved even 
for finite and extremely short GNRs, regardless of the 
level of theory used. 
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